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ABSTRACT 

The initial mass function (IMF) of metal-free stars that form in the initial 
starburst of massive (virial temperatures > lO^K) metal-free protogalaxies is 
studied. In particular, we focus on the effect of II2 photodissociation by pre- 
existing stars on the fragmentation mass scale, presumedly determined by the 
Jeans mass at the end of the initial free-fall phase, i.e., at the so-called "loiter- 
ing phase," characterized by the local temperature minimum. Photodissociation 
diminishes the Jeans mass at the loitering phase, thereby reducing the fragmenta- 
tion mass scale of primordial clouds. Thus, in a given cloud, far ultraviolet (FUV) 
radiation from the first star, which is supposedly very massive (~ lO^M©), re- 
duces the mass scale for subsequent fragmentation. Through a series of similar 
processes the IMF for metal-free stars is established. If FUV radiation exceeds 
a threshold level, the star-forming clumps collapse solely through atomic cool- 
ing. Correspondingly, the fragmentation scale drops discontinuously from a few 
time IOM0 to sub-solar scales. In compact clouds (< 1.6kpc for clouds of gas 
mass 10^ Mq), this level of radiation field is attained, and sub-solar mass stars 
are formed even in a metal-free environment. Consequently, the IMF becomes 
bi-modal, with peaks at a few tenths of Mq and a few times IOMq. The high- 
mass portion of the IMF, ^high("^*)5 is found to be a very steep function of the 
stellar mass m*, ^high("^*) oc m~^. Therefore, the typical mass scale of metal- free 
stars is significantly smaller than that of the very first stars. In the Appendix 
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we study the thermal instabihty in collapsing primordial prestellar clumps, and 
discuss why the thermal instability occuring during the three-body H2 forma- 
tion does not appear to manifest itself in causing further fragmentation of such 
clumps. 

Subject headings: cosmology: theory — galaxies: formation — stars: formation 

1. Introduction 

Recent theoretical studies claim that the first stars are very massive, typically about 
a few times IO^Mq to IO^Mq, on the basis that (i) the fragmetation mass scale owing to 
H2 cooling is large, about IO^Mq (Bromm et al. 1999, 2002; Abel ct al. 2000) and (ii) the 
absence of dust grain formation allows most of the fragment mass to accrete onto a protostar 
without being halted by stellar radiation feedback (Omukai & Palla 2001, 2003). 

The most iron-deficient star yet observed, HE 0107-5240 ([Fe/H]=-5.3, = O.SMq] 
Christlieb et al. 2002), however, presents a puzzle in this regard (see Beers 2003 for a 
concise summary of the pertinent issues). Shigeyama, Tsujimoto, & Yoshii (2003) suggest 
that HE 0107-5240 was most probably born as a metal-free Population HI star, and has been 
metal-enriched by accretion of interstellar gas, rather than being born as a Population H 
star. Whether or not this picture is correct, the very existence of HE 0107-5240 demonstrates 
that subsolar-mass star formation is possible even in metal-free or extremely metal-poor 
environments. If HE 0107-5240 was born as a Population II star, the fragmentation mass 
scale can clearly be reduced to < Mq as long as a significant amount of metals has been 
depleted into dust grains (Schneider et al. 2003). Even so, if a completely metal-free star is 
found in the future, another solution will be required. This possibility has motivated us to 
examine whether it is feasible to form subsolar-mass stars from primordial gas. 

While the first stars were probably very massive, it is important to consider whether 
all metal-free stars must be so massive. Note that metal-free stars might not be exclusively 
associated with the very first stars, although theoretical studies often make this assumption 
(e.g., Nakamura & Umemura 2001; Abel et al. 2002; Bromm et al. 2002). The key is that, 
depending on the nature of Population HI stars, the fragmentation process of primordial 
clouds can be significantly altered owing to stellar feedback, such as the UV radiation field, 
the kinetic effect of stellar winds and supernovae, etc. In fact, in the presence of a strong 
far ultraviolet radiation field, metal-free star-forming clumps can collapse solely by atomic 
cooling, with H2 formation prohibited by photodissociation (Nakada & Yoneyama 1976; 
Hasegawa, Yoshii, & Sabano 1981; Omukai 2001). Then, the fragmentation mass scale can 
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be reduced to sub-solar scales. 

In this paper, we study the evolution of a pregalactic cloud, including its star formation 
history and the growth of a UV radiation field within it. In particular, we focus on conditions 
suitable for low-mass star formation, and consider the initial mass function (IMF) of metal- 
free stars. 

Small halos whose virial temperature is less than 8000K cool only by H2 rotational/ vibrational 
line emission. After the very first episode of star formation, the entire mass of molecular gas 
is photodissociated by a single or a few massive stars, and subsequent star formation activity 
is quenched as long as the photodissociating star is present (Omukai & Nishi 1999; Ciardi, 
Ferrara, Governato, & Jenkins 2000; Glover & Brand 2001; Machacek, Bryan, & Abel 2001). 
In order for star formation to continue, the pregalactic cloud must be massive enough, i.e., 
have virial temperature > 8000K or correspondingly, a total mass > IO^Mq. In such clouds, 
which we consider in this paper, metal-free star clusters will form eventually as a result of 
succesive star formation processes. 

The IMF of those metal-free star clusters is of particular importance to the subsequent 
evolution of the universe, e.g., the metal enrichment, reionization, etc. Consequently, the 
metal-free stellar IMF has been considered by a number of previous authors. Among them, 
Yoshii & Saio (1986) studied it in the framework of Silk (1977) 's IMF theory, which assumes 
the opacity- limited fragmentation scheme of Hoyle (1953), whereby the fragmentation is 
supposed to continue until the clump becomes optically thick. However, according to recent 
numerical simulations, the fragmentation mass scale may instead be determined when the 
isothermality breaks, before the opacity becomes important, i.e., at the epoch when the 
clump temperature begins to increase (Bromm et al. 2002). 

By analysing the fragmentation process of filamentary primordial clumps, Nakamura & 
Umemura (2001) and Bromm et al.(2002) obtained similar estimates for the fragmentation 
epoch of low-density filaments (nn ^ lO^cm"^). Such filaments fragment with masses of a 
few times lOOM©. Higher density filaments (nn ^ lO^^cm^'^) fragment, with masses on the 
order of IMq, when the clump becomes optically thick to H2 lines (see also Uehara et al. 
1996; Uehara & Inutsuka 2000). On these grounds, they claimed that the metal-free IMF 
becomes bi-modal, with peaks around a few times lOOM© and IMq, assuming the initial 
density of filaments are distributed homogeniously up to > lO^cm"^. However, it is not clear 
whether the initial densities extend to such high values in realistic situations. Indeed, the 
3-D numerical simulations by Abel et al. (2000; 2002) and Bromm et al. (1999; 2002) have 
not succeeded in producing such low-mass fragments. 

Above all, we emphasize that much of the recent work on the metal-free IMF does not 
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take proper account of feedback effects from already-formed stars. Although various forms 
of feedback effects can be pointed out, we consider here only the contribution of the FUV 
radiation field arising from pre-existing stars as a first step toward solution of this problem. 

This paper is organized as follows. In §2, the effect of a FUV radiation field on the 
fragmentation mass scale of primordial gas is discussed. In §3, the evolution of a pregalactic 
cloud, and the mass spectnim of formed stars are discussed using the result obtained in §2. 
In §4, a number of important parameters for pregalactic chemical evolution are presented. 
A brief summary and discussion is then presented in §5. In the Appendix, the possibility of 
fragmentation due to thermal instability is re-examined. 

2. Fragmentation of Primordial Gas Clumps in a Far Ultraviolet Radiation 

Field 

In this section, we study the effect of FUV radiation on the fragmentation mass scale of 
primordial gas. We assume that a pregalactic cloud has an inhomogeneous density structure, 
consisting of individual "clumps". As these clumps collapse, they eventually fragment into 
dense cores ("fragments") of mass mfrag- We consider here the evolution of each clump, in 
particular its fragmentation mass scale. 

Although the detailed dynamical evolution of such clumps is rather complex, in general, 
clumps more massive than several times the Jeans mass are thought to first collapse in a 
disk-like manner, and subsequently fragment into filamentary clumps (Miyama, Narita, & 
Hayashi 1987a,b). Finally, with further collapse, the filamentary clumps fragment again into 
dense cores (Inutsuka & Miyama 1997). All these processes proceed approximately on the 
free-fall timescale. The density inhomogeneity is amphfied by self-gravity in the course of the 
collapse. The central region experiences run-away collapse, but leaves the outer envelope of 
lower density unevolved. The size, namely the diameter in the spherical and filamentary case 
or thickness in the disk case, of the central region is about the Jeans length Aj, regardless 
of the shape of the clumps. 

The thermal and chemical evolution of primordial clumps irradiated with a FUV radi- 
ation has been investigated by Omukai (2001). Here, we use the same model as in Omukai 
(2001), with minor modifications so that the collapse timescale is changed to the free-fall 
timescale for the filamentary configuration. Hence, we adopt t^oi = l/\/4Gp (i.e., Inutsuka 
& Miyama 1997) rather than the equivalent expression for the spherical case, 
where p is the central density of clumps. This change has only a small effect on our results. 

The overview of the model is as follows. In this model, the thermal and chemical 
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evolution for the central region of the clumps is computed. Within the central region, whose 
length scale is taken as Aj/2, physical quantities are assumed to be homogeneous. The 
density increases at the free-fall rate: 

dt tcoi 

The thermal evolution is followed by the energy equation: 

dC/th P dp 



dt p2 



~ ^netj (2) 



where Uth is the thermal enegy per unit mass, and P is the thermal pressure. The net cooling 
rate, >Cnet, consists of the contributions from the radiative cooling of atomic hydrogen lines, 
molecular hydrogen lines, continuum radiation from primordial gas, and the cooling/heating 
associated with chemical reactions. The optical depth is evaluated by the length scale of the 
central region Aj/2. In the case that the clump is optically thick, absorption and scattering 
effects are taken into account by the escape probability method. The radiation field from 
outside the clump is also attenuated with the same optical depth. The main effect of radiation 
field is photodissociation of molecular hydrogen, which is the major coolant in primordial 
gas clumps. Non-equilibrium chemistry between H,H2,H+,H^, and is solved. Helium 
is neglected since it is thermally inactive in the considered temperature range. For more 
details, we refer the interested reader to Omukai (2001). 

The radiation field is assumed to be produced only by stars inside the pregalactic cloud, 
while the possible contributions either from the extragalactic background radiation or from 
active galactic nuclei are neglected. As discussed §3 below, only very massive (40 — 2OOOM0) 
stars are considered as sources of this UV radiation. Since the effective temperature of 
such massive stars increases only moderately with stellar mass, for simplicity, we assume 
that the radiation field inside the cloud is approximated by a diluted blackbody at a single 
temperature T^ad and take Trad = lO^K, which is a typical value of metal-free very massive 
stars (Ezer & Cameron 1971). Also, the ionizing photons are assumed to be trapped in the 
Hll region formed around the massive stars, and not to contribute to the average radiation 
field in the protogalaxies. Then, using a parameter W, the mean intensity can be expressed 
as 

J,^WB,{T,^), (3) 

for hu < 13.6eV, and Ji, — for hu > 13.6eV. The parameter W is called the dilution factor, 
hereafter. 



Figure 1 shows the temperature evolution of collapsing primordial clumps in the presence 
of FUV radiation, plotted as a function of the density. The dilution factors, W, are as 
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listed in this figure. The initial temperature is taken arbitrarily to be 300K at O.lcm"^, 
because the evolution at higher density, where the fragmentation occurs (at H number density 
riH ^ lO^cm"^), is not afi^ected by this choice. 

First, we describe the temperature evolution for primordial clumps without and external 
radiation field (e.g., Saslow & Zipoy 1967; Matsuda, Sato, & Takeda 1969; Yoneyama 1972; 
Yoshii & Sabano 1979; Carlberg 1981; Palla, Salpeter, & Stabler 1983; Omukai & Nishi 
1998). As seen in Figure 1 (the curve labelled W = 0), the temperature of the clump first 
rises adiabatically, then starts decreasing abruptly by forming sufficient H2 for cooling. The 
temperature decrease stops at riu ~ 10^~^cm~^, where rotational levels of molecular hydrogen 
reach the local thermodynamic equilibrium. After that, the temperature increases gradually. 
The actual temperature increase in this phase is slower than shown in Figure 1, owing to 
the slower collapse rate until the formation of nearly spherical objects by fragmentation. We 
call this slow collapse epoch the "loitering phase" (e.g., Bromm et al. 2002). 

The loitering phase clearly plays a role in the fragmentation of primordial clumps, and 
has been studied by Bromm et al. (1999, 2002) and Abel et al.(2000, 2002) by way of 3-D 
hydrodynamical simulation, and by Nakamura & Umemura (2001) by 2-D hydrodynamics 
assuming that filamentary configurations are formed. According to these studies, the Jeans 
mass at the loitering phase is imprinted as a typical mass scale for fragmentation. Fragmen- 
tation occurs around this point (rin ~ 10^~^cm~^), and the typical fragment mass scale is 
the Jeans mass then. 

Next, let us see the effect of FUV irradiation on the thermal evolution. Under the 

presence of a UV radiation field, the formation epoch of H2 is delayed until higher densities 
and temperatures are reached (see Figure 1). The density and temperature at the loitering 
phase also increase. The net efi^ect of this is to reduce the Jeans mass at the loitering phase. 
In other words, the stronger the FUV radiation, the smaller the fragmentation m,ass scale. 

If the radiation intensity is even more increased, say to the level of > 1 x 10~^ 
in terms of the dilution factor, the photodissociation (for uh ^ IC^cm"^) and collisional 
dissociation (for higher densities) prohibit formation of sufficient H2 for cooling during the 
entire prestellar collapse phase. These clumps cool only through atomic radiation processes, 
namely, H line cooling for nn < lO^cm"^ and H~ bound- free transitions for higher densities 
(see the evolutionary path for 1^ = 1.1 x 10~^ in Figure 1). The evolutionary paths for the 
clumps with W > 1.1 x 10~^ are almost identical to one another except for the shallow "pit" 
around nn ~ lO^cm"^ for 1^ = 1.1 x 10~^, which dissappears in the presence of stronger 
radiation. An important aspect of the collapse by atomic cooling is that the temperature 
continues to decrease up to rin ~ lO^^cm^'^ Consequently, fragmentation occurs at very 
high density, and sub-solar mass scale fragments are formed. 
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Next let us evaluate the fragmentation mass scale inferred from the physical conditions 
at the loitering phase. We define the loitering point as the instance of the temperature 
minimum around the loitering phase. Fragmentation does not take place immediately at the 
loitering point. Instead, the clump contracts more, by about 2 orders of magnitude in density 
(e.g., Nakamura & Umemura 2001), owing to inertia before the fragmentation. Taking 
this effect into consideration, we evaluate the density and temperature at fragmentation as 
follows. After the loitering point, the clump is assumed to fragment by contracting 2 orders 
of magnitude in density. During this period, the temperature is assumed to be constant for 
molecular cooling clumps, while for atomic cooling clumps the temperature rises (as in Figure 
1) because the latter become optically thick before fragmentation. Using the density and 
temperature at the time of fragmentation, the fragment mass scale is estimated as follows. 
The wavelength of maximum growing mode for fragmentation is Amgr = 27ri?corc) where the 
core radius of the filament -Rcorc = -^1/2. Thus, the typical mass scale for fragments is given 
by mfrag = p7ri?po;,gAmgr, whcrc p is the density at the fragmentation (e.g., Uehara et al. 



In order to collapse and fragment, the mass of the clumps must be more massive than 
the maximum Jeans mass attained during the collapse. Therefore, only sufficiently massive 
clumps can give a birth to dense cores and stars, eventually. Probably, in the pregalactic 
clouds, some clumps satisfy this condition and some do not. Note that, in our model, as 
long as the clumps are massive enough, the fragmentation mass is only a function of the 
radiation field since we have adopted the same thermal evolution given in Figure 1 for all 
those clumps. 

Fragmentation mass scales computed in this way are shown in Figure 2 as a function 
of FUV radiation intensity. In the case of no radiation field {W — 0), the fragmentation 
mass is mfirst — 2OOOM0, which is consistent with the numerical simulations (e.g., Bromm 
et al. 2002). The fragmentation scale decreases gradually from the initial value mgrst as 
the radiation intensity increases. When the radiation exceeds the threshold level, VFnoH2 — 
1.1 X 10^^^, a transition from the Il2-cooling triggered collapse to the atomic-cooling one 
occurs. Correspondingly, the fragmentation scale drops discontinuously from the transiton 
mass, mtr — 4OM0, to sub-solar scales, miow — O.SM©. 

For an insight into possible values for W ^ let us consider here the radiation field at 
distance r from a star of radius i?* and of effective temperature Teflf- The mean intensity 
there is 



1996). 




(4) 
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The dilution factor is thus given by 

Using the relation L = AnRlcrT^Q, and noting the luminosity of massive stars L is close to 
-^Edd; the dilution factor can be expressed as 



where m* is the mass of the star. Then, if the massive stars (> SOMq) are predominantly 
formed and their mean separation is ~ a few pc, the threshold level VFnoH2 — 1-1 x 10~^^ 
can be reached. In §3, we discuss whether this level of radiation field is reached during the 
evolution of pregalactic clouds. 

In this paper we assume that the dense clumps formed by fragmentation around the 
loitering phase do not produce sub-fragments. The possibility of sub-fragmentation due to 
thermal instability has been clamed by some authors (Sabano & Yoshii 1977; Yoshii & Sabano 
1979; Silk 1983; Yoshii & Saio 1986). In the Appendix, we re-examine this question via a 
linear perturbation analysis of the gravitationally collapsing primordial gas. Although the 
thermal instability does occur in the active phase of the three-body H2 formation reaction, 
the instability remains weak and does not appear to cause further fragmentation of the 
clumps. Therefore, we do not consider the sub-fragmentation by the thermal instability in 
the following. 



3. Star Formation in Pregalactic Clouds 
and Establishment of the Stellar Mass Spectrum 

In §2 we discussed the evolution of subclumps in a pregalactic cloud, and obtained the 
fragmentation mass scale as a function of the FUV radiation intensity. Here, let us consider 
the evolution of pregalactic clouds as a whole. The history of the build-up of the stellar 
radiation field and the consequent IMF within a pregalactic cloud is examined. As discussed 
in §2, sub-solar mass stars are formed if the FUV radiation field exceeds the threshold level, 
which is VFnoH2 — 10~^^ in terms of the dilution factor. The condition of the cloud when this 
transition actually occurs is also examined. 
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3.1. A Model for Pregalactic Clouds 

Suppose that pregalactic clouds consist of two-phase medium, i.e., diffuse gas and denser 
clumps. The clumps collapse and form stars eventually, as described in §2. Since the physical 
condition of the clumps, e.g., initial density, temperature, etc., varies from clumps to clumps, 
the star formation epoch in those clumps would differ each other. The star formation rate 
in clumps, local radiation field, and other local quantities, do depend on the local condition 
in the clumps. Here for simplicity, we assume that the quantities averaged over the entire 
cloud, e.g., the averaged star formation rate, are controlled by averaged global quantities of 
the cloud, e.g., the averaged gas density. Specifically, we use a one-zone model of pregalactic 
clouds, where the physical quantities such as the gas density, Pgas, stellar density, p*, and 
stellar radiation field, Wrad, are taken to be averaged values (see e.g., Yoshii & Saio 1985). 

Let the initial gas mass be Md, and the length scale be I. The initial gas density is 
given by 



Pga.(0) = (7) 

The length scale of the cloud / is assumed to be constant in time. While the fragments 
themselves continue their own collapse, the parental pregalactic cloud is kept in dynamical 
equilibrium by the kinetic energy of the fragments. 

The basic assumptions of this model are (i) star formation proceeds on the free-fall 
timescale, and (ii) the stellar mass scale coincides with the fragmentation scale, which is a 
function of the radiation intensity, using the relation in Figure 2. 

The average gas density, pgas, in the cloud is followed by 



dpgas _ Pgas 

dt tsi 

For the star formation timescale, tst, we take the free-fall timescale, tco\ — l/^/iGp^. We 
neglect the recycling of gas ejected by either supernovae or stellar winds, assuming that it 
directly escapes from the cloud. Also, gas is always assumed to be primordial and chemical 
enrichment effecte are not taken into account, since we are now interested only in metal- free 
stars. 

Let the number of stars per unit volume at time t in mass range dm* be ^{'m^,t)dm^. 
By summing up the radiation from all stars, the average radiation enegy density in the cloud 
is given by 
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Urad ^ ~J C{n^*,t)L{m^)dm^. (9) 

The values of luminosity L{m^) of a Population III star with mass m^, are taken from Schaerer 
(2002). Here, the radiation spectrum is supposed to be a diluted blackbody of Trad — lO^K, 
since we consider the radiation from massive metal-free stars (> rritr — 40Mq). The mass 
scale of fragmentation at time t, mfi.ag(t), is determined by Mrad(^) by the relation shown in 
Figure 2. Stars are assumed to turn off after their main sequence lifetimes. Here, the main 
sequence lifetime of stars without stellar wind given by Schaerer (2002) is used. For stellar 
masses higher than 500Mq, whose lifetime is not given in the literature, the same lifetime 
as that for SOOMq is used. 

With the above assumptions, the equation for stellar number density, ^(m*,t), is for- 
mally written as 



dt 



5{m^ - mf^a.g) - ^{rrid, ta) 



dt 



rrid), 



(10) 



where ma is the mass of stars that dies at time t, and td is their formation time, hence 
td = t — t^,{md)■ In the equation above, the first term on the right hand side is the star 
formation rate, and the second term is the stellar death rate. Equations (8) and (10) are 
integrated until Pgas falls below 10 % of the initial value, Pgas(O), or the FUV radiation exceeds 
the threshold value for atomic-coohng collapse, for which we use the threshold dilution factor 
W^noHa = 1 X lO""*^^. Note again that this is still a quahtative, rather than a quantitative, 
approach to study early star formation in pregalactic clouds, and it is of course not realistic 
that only 10 % of the initial gas remains into the ISM. The evolution of pregalactic clouds 
is not affected by this value: only the termination epoch of the calculation is changed. 



3.2. The Conditions for Low-Mass Star Formation 

Two cases of cloud mass Md are studied: (a) = lO^M©, and (b) IO^Mq. For case 
(a), the length scale, I, is varied from 100 to lOOOpc, while for case (b), / = 300 — 3000 pc. 
The evolution of the average radiation field,Mrad) is presented in Figures 3 (for case a) and 4 
(for case b). 

We see in Figures 3 and 4 that the radiation field in relatively compact clouds, whose 
length scale / < 500pc for — IO^Mq or lOOOpc for = lO^M©, reaches the critical level 
for low-mass star formation, while it does not in more diffuse ones. In pregalactic clouds 
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where low-mass stars are formed, the value of W is expected to remain at the threshold level 
WnoH2 until early-formed massive stars begin to die about t ~ 2 x lO^yrs, since low-mass 
stars do not emit FUV radiation. On the other hand, in larger clouds the radiation enegy 
density -Urad fails to reach the critical level 'UnoH2 = ^noH2o7"4d) ^.nd declines after the first 
stars begin to turn off. In other words, there is a threshold cloud length, /tiu below which the 
radiation intensity reaches the critical value. This threshold length is larger for the cloud of 
IO^Mq (/th > lOOOpc) than for that of IO^Mq (/th < lOOOpc). 

We now discuss analytically the behavior of the radiation field and the threshold cloud 
length. For simplicity, we suppose that the gas depletion timescale, tjcp — 4f(0), is longer 
than the Ufetime of the first star, t*(mfirst)- A similar discussion holds in the opposite case 
with some modification. Also, since the stellar lifetime and luminosity-to-mass ratio is only 
weakly dependent on mass for massive stars, in the following analytical discussion we use a 
constant value of t^^high — 2 x 10^ yr and (L/M)^ = 3 x 10'^{L/M)q, which are the values 
for a star of IOOOMq, for all massive stars (> mtr). 

Before significant gas depletion takes place, ai t < t^ep: the average gas density Pgas — 
Pgas(O) and the star formation rate is approximately constant at Pgas(0)/tsf(0). The stellar 
mass density is 



p (t) ~ I Pgas(0)V^sf(0) t < t*,high, ^^^^ 
i Pgas(0)i*,high/^sf(0) i*,high < t < tdep- 

After the gas is depleted significantly, the star formation rate drops to Pgas(^)/4f(^)- Thus 



p*{t) ~ Pgas(^)^*,high/^sf(^) t > Uep- (12) 

By the way, the radiation density is given by 

Mrad(^) = -p*{t){L/m)^, (13) 
c 

as long as all the stars are massive, namely Wrad < MnoH2- Since {L/m)^. is only weakly 
dependent on m*, the radiation field is proportional to the stellar mass density p*. Therefore, 
the radiation intensity first increases linearly in time {t < i*,high), then remains constant 
(^*,high <t < ^dep), and finally declines as the gas in the cloud depletes {t > t^ep)- 

The radiation dilution factor W, which is related to the radiation energy density Urad 
by 1^ = Urad/7.56 X 10^(erg/cm^), is 
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and reaches maximum at t 
than the threshold length 

Zth = 1.6kpc(Mel/10«Mo)=^/^ (15) 

the dilution factor W reaches VFnoH2 — 1 x 10^^^, and sub-solar mass stars are formed 
thereafter. This value of /th agrees well with the numerical result presented in Figures 3 and 
4. 

The number of low-mass stars formed can be estimated by the following considerations. 
Since low-mass stars do not emit UV radiation, u^ad remains constant after reaching ttnoH2 
at t = timsf; where 



iiki) iwk) (j^) '^"^ 

= i*,high- Therefore, only in clouds whose length scale is smaller 



iimsf = 4 X lOV (^^] , (16) 

until the massive stars start to die at t = t*,high- Then, the duration of low-mass star 
formation is t*,high — ^imsf- Multiplying by the star formation rate, the mass density of low- 
mass stars is 



P*,low — , /^N (^*,high ^Imsf)- (l*^) 

Metal-free star formation continues for the duration At = t*,highi i-e., from the first star 
formation until the massive stars begin to explode as supernovae. At time t = timsf in 
this period, low-mass star formation starts. Since the star formation rate is approximately 
constant at Pgas(0)/tsf(0), the mass fraction of metal-free stars that is locked into low-mass 
stars is given by the ratio of the two durations: 



flovf 



P*,low ^*,high ^Imsf 

P*,Z=0 ^*,high 

l-0 2(^V7^^"^'^ 

' IkpCy/ VlO^^O 



1 - 



(18) 
(19) 
(20) 
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where p*,z=o is the density of metal-free stars. 



3.3. 



The Stellatr Initial Mass Function 



Because of the mass gap between low-mass stars (~ Q.SMq) and high-mass stars (> 
4OM0), the metal- free IMF becomes bi- modal. 

In our idealized model, the low-mass IMF is delta-function-like: only stars of miow — 
0.3Mq are formed (see Figure 2). Actually, a distribution of the fragmentation mass around 
this value would be present. Competitive accretion of ambient gas would create some distri- 
bution of the protostellar mass. Also, the coagulation between fragments might form more 
massive cores. As a result of these processes, the low-mass IMF would be formed, although 
we do not discuss it further here. Recall that the peak ratio between the high and low-mass 
IMF can be obtained by the mass fraction of low- mass stars /low given by equation (20). 

On the other hand, the high-mass IMF, ^high(?7i*), of metal-free stars that are more 
massive than the transition mass mtr — 4:0Mq has been obtained in our model, and is shown 
in Figure 5 for clouds of Md = lO^M© (case a). In three cases of I — 300, 500, and lOOOpc, 
the shape of the IMF is identical, and the value of IMF at the same mass is proportional to 
l-^. In the case of / =3 kpc, the radiation field does not reach the threshold value, and no 
low-mass stars are formed (see Figure 4). In this case, the instantaneous mass of forming 
stars first reaches a minimum of I4OM0 and then increases as the radiation intensity declines. 
The stars of < 3OOM0 can be produced at two different epochs, first when the radiation 
density is increasing and second when the radiation density is decreasing. The stars formed 
at different epochs are shown separately in Figure 5. The actual IMF is the sum of them. 

We now present an analytic derivation of the high-mass IMF, ^high, following Silk (1977). 
Here, we consider stars formed at t < t*(mfirst), namely, during the time when the first star 
is still alive, as those formed later are probably metal-polluted. Thus, the maximum mass of 
existing stars is just the mass of the first star, mfirst- The radiation energy density at time t 
is given by 




(21) 



where m*{t) is the mass of forming star at t. Differentiating both sides of the equation (21) 
with m{t), we obtain the mass spectrum 
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C (iWrad 



(22) 



lL{m^) dm* 



Suppose now the relations i^rad and L oc m^. In this case, the IMF becomes 



^high oc 



a-/3-l 



(23) 



This relation was first derived by Silk (1977). From Figure 2, the approximate relation 
between fragmentation mass and radiation density, mfrag oc W~^^^, holds for massive stars 
of m^: ~ 1OOM0. For these massive stars, the luminosity-mass relation L oc holds. Then, 
it follows that a = —3 and /? = 1 for those stars. Substituting these values into relation 
(23), we finally obtain the metal-free IMF for massive stars: 



The mass spectrum becomes a very steep function of mass because of the negative value 
of a. The normahzation constant for the relation (24) depends on the length scale of the 
pregalactic cloud as oc as seen in equation (21). Thus, using the values in Figure 5, the 
IMF can be expressed apprximately as 



Note that the IMF is independent of the mass of the pregalactic cloud (see eq. 21). Also the 
transition mass mtr — 4OM0 does not depend on the cloud mass because this comes from 
the discontinuity in the mfrag — W relation (Figure 2), which is independent of the cloud 
mass. 

A bi-modal IMF for metal-free stars has also been proposed also by Nakamura & 
Umemura (2001), as mentioned in §1. It is interesing to note here the differences between 
our IMF and theirs. They pointed out that the IMF becomes bi-modal, with peaks around 
1 and 100 Mq, if dense (nn > lO^cm"'^) filamentary clumps formed in the course of collapse. 
It is, however, still uncertain that such dense filaments are formed in reality. On the other 
hand, we consider here only the low-density filaments that fragment initially during the loi- 
tering phase. Even so, the sub-solar mass peak appears, owing to photodissociation feedback 
from pre-existing stars. Recall also that the production mechanism of the low-mass stars is 
different in the two treatments. Nakamura and Umemura (2001) considered fragmentation 
of dense molecular-cooling clouds around ~ lO^^cm"^, while we considered fragmentation 
of atomic-cooling clouds around ~ lO^^cm"^. 



Chigh oc m/. 



(24) 




(25) 
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4. Metal Enrichment by First-Generation Stars 

The chemical enrichment by the first-generation stars has particular importance for 
the subsequent chemical evolution of protogalaxies. Also, the masses of those stars can be 
constrained by observing the abundance ratios of elements in second-generation stars. In 
§3, the IMF for metal- free stars was predicted to be bi- modal with peaks around O.SMq and 
AOMq. In this section, we discuss the chemical yield of the first-generation stars under the 
predicted IMF. 

Since low mass stars do not eject metals (until after completing their long lifetimes), we 
do not need to know the detailed shape of the low-mass IMF. The mass ratio of low-mass 
stars, /low, given by equation (20), suffices for our purpose. For simplicity, we approximate 
the high-mass portion of IMF as 

Chigh oc (mtr < < murst) (26) 

in this section. We define the mean mass fraction of ejected element X by super novae 

(27) 

and the mean remnant mass fraction 

«high = J C?"7,*mremChigh/ J c?m*m*^high- (28) 

Using these quantities, the yield of element X by the first-generation stars is expressed as 

px = ^^'h'gh(-^ ~ -^1°^) (29) 

/low + (1 ~ /low)Q;high 

Considering the uncertainty of the transition mass, mtr, we calculate the yield for two 
cases of mtr = 20 and AOMq, respectively. We fix mgrst — 2OOOM0 since the uncertainty of 
this causes httle difference owing to the steepness of IMF. The values of ex,high and ahigh are 
presented in Tables 1 and 2 for mtr = 20 and AOMq, respectively. We use the Population III 
supernova (SN) yields by Umeda & Nomoto (2002). Metals are ejected by stars in the mass 
range IS-SOMq by Type II SNe, and 150-270Mq by pair-instability SNe. Metal-free stars 
whose mass is between these ranges, as well as those more massive than 270 Mq, collapse into 
black holes with little metal ejection. Since the fraction of gas returning to the interstellar 
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medium is quite uncertain, we assume here that all of the stellar material is eventually locked 
into black holes in this mass range. 

A comparison of predicted elemental abundance ratios with the observed most iron- 
deficient star HE 0107-5240 (Christheb et al. 2002) is shown in Table 3. When the cutoff 
mass, mtr, is higher than SOMq, metal enrichment proceeds only via pair-instability SNe, and 
the abundance ratios become identical to this single source. Christlicb ct al. (2002) claimed 
that a Type II SN of a 20 — 25Mq metal-free star reproduces the observed abundance pattern 
heavier than Mg, while lighter elements could be accounted for by self-enrichment or mass 
transfer from a more massive AGB companion star. On the other hand, Schneider et al. 
(2003) pointed out that the abundance pattern from a pair-instability SN of a 200 — 22OM0 
star also agrees with the limited number of elemental abundances presently available for 
this star. So far, the enrichment either by Type II or by pair-instability SNe seems to 
be consistent with the observed abundance ratios. As a result, the abundance pattern of 
HE 0107-5240 heavier than Mg is well-reproduced in all cases shown in Table 3 within 0.2-0.3 
dex, although in the case of mtr > 3OM0 the amount of Mg relative to Fe might be too small. 

Recently, Umeda & Nomoto (2003) proposed the "mixing and fallback" scenario: the 
processed material is assumed to be mixed uniformly in the region from I.SMq to 6.OM0, 
and almost all material in this layer falls back to the central remnant; only a small fraction is 
ejected from this region. Although this beautifully explains the abundance ratios of HE 0107- 
5240, inchiding the large C and O abundances relative to Fe, this theory contains unknown 
free parameters, namely the extent of mixing and the mass-cut. 

Note, however, that we do not take it as yet established that HE 0107-5240 is a second- 
generation star formed from metal-enriched material. We stress here that since sub-solar 
mass stars can form in a metal-free environment, HE 0107-5240 might be a first-generation 
star metal-polluted due to subsequent accretion (Yoshii 1981; Shigeyama et al. 2003). The 
large [C/Fe](=+4.0) ratio of HE 0107-5240, as well as the large [0/Fe] (=+2.4) ratio (private 
communication; M. S. Bessell), can be explained by the self-enrichment in the AGB phase 
(private communication M. Y. Fujimoto; Fujimoto, Ikeda, & Iben 2000). The accretion 
scenario demonstrates that the gap in the metallicity distribution function of observed low- 
metallicity stars, between [Fe/H]=— 5.3 and the current lower limit [Fe/H]= -3.5 to -4.0 of 
extreme metal-poor stars could be real, as shown in Figure 2 of Shigeyama st al. (2003). 

On the other hand, if the explosion of massive metal-free stars was responsible for the 
birth of the [Fe/H]=-5.3 star as a second-generation star (Schneider et al. 2003; Umeda 
& Nomoto 2003; Bonifacio, Limongi, & Chieffi 2003), there must exist many more second- 
generation stars having a wide variety of metal abundances below [Fe/H]=-3.5. Therefore, 
this massive Population HI scenario prefers the existence of a metallicity distribution without 
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any appreciable gap below [Fe/H]=-3.5, and future large samples of extremely metal-poor 
stars will clearly discriminate between the competing scenarios proposed thus far. 

5. Summary and Discussion 

We have studied the initial mass function (IMF) of the first-generation stars under the 
assumption that the fragmentation of prestellar clumps occurs around the "loitering phase," 
i.e., when the temperature stops decreasing. The fragmentation mass scale decreases as the 
far ultraviolet radiation field becomes intense. If the very first stars are very massive and 
create intense radiation fields inside their parent clouds, stars formed later are smaller than 
those formed earlier. The stellar initial mass function is established by this mechanism. 
Low-mass star formation occurs only in clouds whose length scale / is less than the threshold 
length /th ^ 1.6kpc(Mci/lO^M0)3/7, where Md is the initial gas mass of the cloud. In such 
clouds, the IMF becomes bi-modal: the high-mass peak is about 40Mq and the low-mass one 
locates at about 0.3Mq. Numerous low- mass metal-free stars are formed, greatly exceeding 
the number of massive meta-free stars. The high-mass portion of the IMF is found to be 
Chigh = 3 X 10~^(stars Mq^ pc~^)(Z/lkpc)~^(m*/5OM0)~^, which is a very steep function 
of stellar mass m*, with a power-law index of —5, in comparison with the Salpeter IMF's 
index of -2.35. Thus, although the very first stars are possibly very massive, ~ IOOOMq, 
the typical mass scale for metal-free stars is smaller, and depends on the length scale of the 
pregalactic clouds. 

If all metal-free stars were indeed as massive as the very first stars, they would evolve 
to black holes without metal ejection. The surrounding gas remains metal-free, and the 
very massive star formation mode continues indefinitely. This is the so-called star formation 
conundrum pointed out by Schneider et al. (2002). Our scenario, where the UV radiation 
from the preceding stars allows production of low-mass stars from metal-free gas, might 
provide a solution to this conundrum. 

Recently, Mackey, Bromm, & Hernquist (2003) studied the star formation history in the 
high-redshift universe, in particular the transition from Population III to Population II star 
formation. They assumed that the same top-heavy IMF is maintained not only in low-mass 
halos that cool only via molecular coohng, but also in halos that are massive enough to cool 
via atomic hydrogen lines. According to our results, however, the mass scale of stars that 
form in massive halos is smaller than those that form in low-mass halos. 

We now understand that low-mass star formation occurs in a primordial cloud only 
if the radius of the cloud is below the threshold kh ^ 1.6kpc(Mci/lO^M0)^/^. We should 
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ask whether this is the case in reahsty. The virial radius of a pregalactic cloud, Tyir — 
2kpc(Mci/lO^M0)^/^(l + ^vir/20)~^, that cohapscs at redshift z^i^. should be regarded as 
the maximum radius of the cloud, since it cools and contracts after virialization to reach 
dynamical equilibrium as a whole. Comparing /th and ryij., we conclude that pregalactic 
clouds with / < /th are not rare objects. Therefore, sub-solar mass metal-free stars are 
generically formed in the framework of stardard cosmology. If metal-free low-mass stars are 
not be found in the Galactic halo, an exotic scenario for structure formation must be invoked. 

Our model for pregalactic clouds is obviously oversimplified, and should therefore be 
improved upon before the above conclusion is taken for granted. In the following, wc sum- 
marize our assumptions and their limitations. More elaborate works taking these points into 
account are desirable in future studies. 

• We have used a one-zone model for pregalactic clouds. Inhomogeneity inside a pre- 
galactic cloud is not taken into account. In particular, we have assumed that the star 
formation timescale is given by the averaged free-fall time. Of course, the evolution of 
pregalactic clouds and initial mass function of first-generation stars crucially depend 
on this assumption. Although this assumption has not been fully justified, there are 
some observational suggestions that the global star formation rate (SFR) in galaxies 
is indeed determined by galactic scale conditions (e.g., Kcnnicutt 1998). 

First, we consider the effect of the spatial ffuctuations of SFR within the pregalactic 
cloud. If the star formation is very active in a region in comparison with the rest of the 
cloud, this star-forming region itself can be regarded as an almost independent cloud. 
Since the length scale of the star-forming region is smaller than that of the entire cloud, 
this allows easier build-up of local radiation field and easier low-mass star fomation. 

Next, we discuss the effects of different values of the global SFR. Possible changes in the 
averaged SFR do not affect the high-mass IMF directly, since the equation (22) does 
not contain the SFR. On the other hand, the condition for low-mass star formation 
discussed in §3.2 will be altered. If the star formation process is slower than that 
assumed here, the radiation field will increase more slowly. This results in a smaller 
value of the threshold length Zth for low-mass star formation than in equation (15). 

• We have only considered the averaged radiation field and have neglected its local 
fiuctuations. In reality, the radiation field is stronger in a region closer to a massive 
star than in outer regions. Therefore, dispersion around the average fragmentation 
scale must exist at any time. This effect might affect the IMF in the pregalactic 
clouds. 
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• We have derived the fragmentation mass scale from the condition at the local temper- 
ature minimum, or "loitering point". Because of this simplified treatment, the value of 
fragmetation mass scale should be regarded as precise only in the order-of-magnitude 
sense. More accurate dctermintation of fragmentation mass scale under UV radiation 
field is an interesting topic for future studies (see e.g., Bromm & Loeb 2003; Yoshida et 
al. 2003). Also, the fragmentation epoch depends on the initial density and line-mass 
of filaments (Nakamura & Umemura 2001). The dense and large hne-mass filements 
collapse until very high density (lO^^cm"^) without fragmenting. It is important to 
specify these properties of filaments in the realistic cosmological context. 

• We have assumed that the stellar mass scale coincides with the fragmentation scale. 
If the accretion rate drops suddenly for massive protostars, as suggested by Abel et 
al. (2002), the stellar mass scale falls short of the fragmentation scale. This depends 
crucially on the late-time behavior of the accretion rate, which has not yet been fully 
explored. 

• We have assumed that the ionizing photons are trapped in the Hll region formed 
around the massive stars, and do not contribute to the average radiation field in the 
protogalaxies. This assumption depends on the size and distribution of the HII regions. 
Prom comparison of the results by Susa & Kitayama (2000) and Omukai (2001), we 
expect that the thermal evolution does not deviate so much from those we used in 
§2 even if the ionizing photons are included. Therefore, even if this assumption is 
invalidated, the results will not be altered significantly. 

• We have neglected kinematical effects of the HII regions formed around massive stars. 
If either the density around the star or the protostellar accretion rate is sufficiently 
high, the HII region is confined to a small region around the exciting star. If not, 
the Hll region expands due to high internal pressure, and the entire pregalactic cloud 
might be disrupted by photoevaporation. If the pregalactic clouds are disk-like, an Hll 
region extends preferentially in the vertical direction to the disk, and it ends up with a 
hole like a chimney (Franco, Tenorio-Tagle, & Bodenheimer 1990). Therefore, even if 
photoevaporation by massive stars occurs, it may not be sufficiently violent to disrupt 
the entire pregalactic cloud. However, more thorough studies are necessary to draw 
this conclusion with confidence. 

• We have assumed that formed stars start shining as main sequence stars immediately. 
Stars are, however, formed by protostellar accretion, and the accretion phase precedes 
the main-sequence phase. If the mass accretion rate is mace, the duration of accretion 
phase for a star of mass m* is iacc = m*/macc- In the case of primordial star formation, 
the accretion rate is very high: 10""^ — lO^^Mg/yr (Stabler, Palla, & Salpeter 1986; 
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Omukai & Nishi 1998). Then, the accretion phase is short in comparison with the 
stellar lifetime. However, for the most massive stars (< IOOOMq), the accretion phase 
lasts for a significant part of the expected stellar lifetime. For such high accretion rate, 
the accretion flow becomes optically thick to H~ bound-free absorption. The effective 
temperature, Tcs, of accreting protostars is about 6000K, while at the main-sequence 
stage it becomes T^s ^ lO^K (Stahler et al. 1986; Omukai & PaUa 2001, 2003). Thus, 
if the accretion phase is prolonged, the emitted FUV radiation decreases significantly. 

• We have assumed that the star formation timcscale is similar to the free-fall time of 
the cloud. In other words, we considered only the spontaneous mode of star formation. 
Triggered star formation by stellar winds or SNe (Elmcgrccn & Lada 1977) has not 
been included in our model, since both stellar winds or pulsational mass losses are 
expected to be rather weak for metal-free stars (Kudritzki 2000; Baraffe, Heger, & 
Woosley 2001). 
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A. Thermal Instability in Collapsing Metal-Free Clumps 

In this appendix, we study the thermal-chemical instability using a linear perturbation 
theory, following Sabano & Yoshii (1977) and Yoshii & Sabano (1979). The energy equation 
(Al) and the equation of H2 formation/dissociation (A6) are perturbed under the isobaric 
condition 5P — 0. The energy equation is 

dU _ Pd^ap ^ _^ 

dt p dt ^ ' 

where U is the internal enegy per unit mass, including the thermal and chemical binding 
energy 

2 iijn^ 2 rriH 

and C — is the cooling rate by H2 line emission. In equation (A2), / is the molecular 
fraction ^, and x is the chemical binding energy of H2. Note that, unlike in §2, the cooling 
rate C does not include chemical cooling/heating. Instead, the definition of the internal 
energy U is modified to include the chemical binding energy. The equation of state for an 
ideal gas is used: 

P = (A3) 
/^-^ = 1 - (A4) 

7h2 - 1 3 6 

where the adiabatic index for molecular hydrogen, = 7/5, has been used, and the electron 
fraction, Xg, has been assumed to be small, Xf. <^1. 



In the equations above. 



and 



The equation for H2 formation/dissociation is 

dt 



f = ^, (A6) 



where J- is the net formation rate for H2. 
By introducing new variables 

X = //-^ = 1 - (A7) 



^/ = 1 corresponds to fully molecular. 
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2mnC 1 1 Y 

and 

e = ^. (A9) 



the equations (Al) and (A6) can be written as 

dlnT 2x (Rnp 
dt x + 2 dt 

and 



+ V = Q (AlO) 



dt 

For perturbed quantities, we use the following notation: 

T p X 



+ Q = 0. (All) 



5t = — , 5p = -, 5x = -- (A12) 



The linear analysis is carried out under the isobaric condition 5P — Q: 

5t + 5p + 5^^ 0. (A13) 

Perturbation of equation (AlO) reads: 

'^^^ (^) II + ^^^^^ + P^''^'' + ^"^^^ - ^) = 0, (A14) 



At 

^ = (A15) 



\x + 2j dt 

where 

~ (x + 2)2''dyn 

We define the dynamical timescale as 

dlnp\ ^ 

bg 

for the collapsing background of the perturbation. 

Perturbation of H2 formation/dissociation equation reads: 

dS 

+ TQt5t + pQp5p + xQA = 0. (A17) 

We assume that 5T,5p, and 5x oc exp(cc;t). Seeking non-trivial solutions of equations 
(A13), (A14) and (A17) for 5t,5p and 5x, we obtain the following equation for the growth 
rate of instabiliby, uj: 

aQUj'^ + aiUJ + a2 = 0. (A18) 
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The coefficients are 



ao = 1 + 



3 iJ, 



ai 



{TCr-pCp-C) 



2kBT 



(A19) 

(A20) 
(A21) 

(A22) 



and 



02 



x + 2 \ 



[{TVt - pVp){xQ, - pQp) - {TQt - pQp) {xV, - pV, - TZ)] (A23) 

(A24) 



prriii 



6 

X 



qiTTr-pJ'p) 



2 ^ A;bT 



{TTt - pT,) 



3t 



(A25) 

(A26) 
(A27) 



dyn 



In the bottom panel of Figure 6, we plot the e-folding growth time of the instabiUty 
te — oj"'^ for a collapsing primordial clump. The collapse rate is assumed to be free-fall, 
^dyn — ts — \J?)'Kl?)lGp. The H2 line cooling rate is computed as in Omukai (2001), where 
the escape probabihty method is used for the coohng by optically thick fines. Also sfiown 
are the background temperature {top) and molecular fraction (middle) for the perturbation. 
We see from this figure that fluctuations grow by thermal instability, in other words, the 
thermal instablity timescale, t^, falls below the collapse timescale for the clump, tg, during 
the active phase of three-body H2 formation process. However, the ratio of the growth rate 
relative to the to collapse rate, ts/te, remains below 2. Since the scale for thermal instability 
is Ath ~ Cgte while the core size is approximately the Jeans length, Aj ~ Cgtff, this means 
that the growing scale is larger than a half of the core size. Therefore, although the thermal 
instability does occur, it may only accelerate the collapse of the central region rather than 
driving further fragmentation. 

Abel et al. (2000) investigated the collapse of a prestellar core by 3-D hydrodynamical 
simulation. In their calculation, the sub-fragmentation during the three-body H2 formation 
has not been observed. They interpreted this as a result of hypothetical turbulent mixing. 
Our argument here may give an alternative explanation of their result. 
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Figures 




Fig. 1. — The temperature evolution for contracting metal- free prestellar clumps as a func- 
tion of FUV radiation density. FUV radiation density is parameterized by the dilution factor 

W — Mrad/o-T^d, where Trad = lO^K. The values of W are denoted in the figure. The collapse 
proceeds by H2 cooling for FUV radiation below the threshold VFnoH2 — 1-1 x 10^^^. Whereas 
for W > VFnoH2) H2 formation is prohibited by photodissocitation, and the clump collapses 
by atomic cooling. The dashed lines indicate the constant Jeans mass. 
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Fig. 2. — The fragmentation mass scale as a function of the dilution factor W, which 
parametrizes the FUV radiation density. 
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Fig. 3. — The evolution of FUV ratiation density as a function of the cloud length scale. 
The initial gas mass of the clouds is — IO^Mq. 



-30- 




Fig. 4. — The same as Fig.3, but for the initial gas mass Md = 10^ Mq. 
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Fig. 5. — The initial mass function of metal-free stars in clouds of Md = IO^Mq. Shown 
are the results for the cloud length scale of 300pc (solid line), 500pc (dotted line), lOOOpc 
(short-dashed line), and SOOOpc (long-dashed line). 
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Fig. 6. — The growth rate of thermal instabihty {t~^ = uj) relative to the free-fall rate along 
the evolutionary path of the collapsing clump. Also shown are the gas temperature and H2 
fraction evolutions used as the unperturbed state of such clump. 
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Table 1: 
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ejected mass 


fract 


ion of element X, 




in the 


case 


of the transition 



mass mtr = 2OM0. The mean remnant mass is ahigh = 0.379 
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Table 2: The same as Table 1, but for the transition mass mtr = 4OM0. The mean remnant 
mass fraction a^igh = 0.984 
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Table 3: Comparison of abundance ratios [X/Fe] with HE 0107-5240 



